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Abstract 

It is assumed that strong chromofields are generated at early stages of ultrarelativistic heavy- 
ion collisions which give rise to a collective deceleration of net baryons from colliding nuclei. We 
have solved classical equations of motion for baryonic slabs under the action of a time-dependent 
longitudinal chromoelectric field. It is demonstrated that the slab final rapidities are rather sensitive 
to the strength and decay time of the chromofield as well as to the back reaction of the produced 
partonic plasma. The net-baryon rapidity loss (5y) = 2.0, found for most central Au-Au collisions 
at RHIC, can be explained by the action of chromofields with the initial energy density of about 
50 GeV/fm 3 . Predictions for the baryon stopping at the LHC are made. 
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It is expected that strong chromofields can develop at early stages of ultrarelativistic 
heavy-ion collisions. There exist different suggestions concerning the space-time structure of 
;hese fields, ranging from string-like configurations as in the color Flux Tube Model (FTM) 
l| to stochastic configurations associated with the Color Glass Condensate (CGC) 2\. This 
picture is most conveniently presented in the c. m. frame where two Lorentz contracted 
nuclei look as thin sheets. After their intersection these sheets acquire stochastic color 
charges as a result of multiple soft gluon exchange. Then strong chromofields are generated 
in the space between the receding sheets. At later times these fields decay into quarks and 
gluons which after equilibration form a quark-gluon plasma. This process has been studied 
by several authors under different assumptions about the field decay mechanism (see e. g. 
Refs. 0,0,3). 

Most previous calculations assume that after interaction the nuclear debris follow the 
light-cone trajectories, thus disregarding their energy losses to produce the chromofield. This 
assumption can be possibly justified only at asymptotically high energies. Furthermore, this 
assumption becomes irrelevant when studying the baryon stopping. Obviously, the energy 
of produced fields and particles is taken entirely from the kinetic energy of the colliding 
nuclei. 

As measured by the BRAHMS collaboration |6], in central Au+Au collisions at highest 
RHIC energy v /J/v r /v = 200 GeV per iVW-pair the baryon energy losses are very significant, 
about 70% of the initial energy. The net-baryon rapidity distributions are substantially 
shifted toward the center of mass from the initial rapidities ±?/o. 

The problem of baryon stopping at RHIC has been addressed recently by several authors. 
In particular, net-baryon rapidity distributions were calculated within the microscopic string- 
based models like UrQMD [7] and QGSM [8]. Although these models implement energy and 
momentum conservation, and thus predict a certain baryon stopping, they are formulated 
in momentum space and do not give a space-time picture of this process. Also, they are 
dealing with hadronic secondaries and therefore preclude the quark-gluon plasma formation. 
Recently, the calculations have been also done within a parton cascade model The 
distribution of valence quarks at central rapidities was also studied in ref. [lo| within a 
QCD motivated approach which however can not be extended to the fragmentation regions. 

In ref. PJJ a simple space-time model was proposed where the baryon stopping was directly 
linked to the formation of strong chromofields. There nuclear trajectories were calculated 
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under assumption that the field is neutralized at a sharp proper time 1 fm/c (see also refs. 
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13]). In this paper we further develop this model and apply it for a more realistic decay 
pattern of the chromofield. Our main goal is to calculate the net-baryon rapidity distribution 
and compare it with the BRAHMS data. 

Following ref. [ijj] we decompose the collision of two ultrarelativistic nuclei into a set of 
pairwise collisions of elementary slabs. At given impact parameter b the positions of the 
target and projectile slabs in the transverse plane are determined by vectors s and b — s, 
respectively. The transverse cross section area of individual slabs is assumed to coincide with 
the nucleon-nucleon inelastic cross section (Inn- The energy and momentum of the projectile 
(a = p) and target (a = t) slabs are parameterized as E a = M a cosh Y a and P a = M a sinh Y a , 
where M a = m_|_iV a is the average transverse mass and N a is the average baryon number of 
slab a. The average baryon transverse mass m± differs from the free nucleon mass mm due 
to internal excitation of slabs acquired at the interpenetration stage. It is parameterized 



in terms of the mean transverse momentum as m± = \Jm 2 N + (p±) 2 . The calculations 
below are made in light-cone coordinates, proper time r = \Jt 2 — z 2 and space-time rapidity 
v = ±ln(*±*). 

nn 

In the Glauber model (see e.g. refs. [15|,116j)the average number of participating nucleons 
from nucleus a in a slab of transverse area <jnn located at a radius- vector s is given by the 
expression 

N a {b, s) = a NN A a T a (b- s) (l - [1 - cr NN Ta(s)} A ^ , (1) 

where a = t for a = p and vice versa, A a is the mass number of nucleus a. The normalized 
functions T a (a = p, t) describe the transverse profiles of the baryon number distribution in 
colliding nuclei. They are obtained by the integration of the corresponding baryon densi- 
ties along the beam direction, A a T a = j p a (f)dz. The average number of nucleon-nucleon 
collisions in an inelastic interaction of two slabs at a radius- vector s is expressed as 

AU(&, s) = a 2 NN A p Tp(b — s)A t T t (s) . (2) 

We assume that the initial energy density stored in the chromofield after collision of two 
baryonic slabs can be parameterized as 

/ s \ a r -"10 
e/(r ;6, s) = e ( — J [A^ co ^(&, s)J , (3) 



where eo is the fitting parameter, second factor describes the collision energy dependence 
with a ~ 0.3, as motivated by small x behavior of the gluon structure function [17|. The 
last factor describes effects of the collision geometry, where exponent (3 is related to the 



apping strings f3 ~ 1.0, but in the case 
18j. In this paper we consider only the 



spatial distribution of the chromofield. For non-over 
of string clustering (3 is expected to be closer to 0.5 
case (3 = 1. 

The equations of motion for a baryonic slab under the action of the longitudinal uniform 
chromofield are obtained by applying the energy-momentum conservation laws across the 
slab. It is assumed that the space between the receding slabs is occupied by the chromofield 
and a partonic plasma, produced by the partial field decay. They exert a certain force on 
the slab from inside. But from the other side the slab has nothing but the physical vacuum. 
This results in a net force acting on the slab. The final system of differential equations 
governing the slab rapidities Y a (r) has the following form (detailed calculations see in ref. 



dM 2 t , . ~ 

T2A(r)P a , (5) 



dr 



= T J^_ (q) 



where P a = M a sinh (Y a - rj a ), E a = ^ M 2 + P 2 , B(r) = a NN (e mc + e f -p) and A(r) = 
°~nn (e P +p)- The plus and minus signs in the r. h. s. of these equations correspond to 
the projectile (a = p) and target (a = t) slabs, respectively. We adopt the initial conditions 
for the slab trajectories ?7a(To) = Y a (r ) = ±y at r =0.01 fm, where ±yo are the initial 
cm. rapidities of colliding nuclei. In expressions above e p and p = c 2 e p are, respectively, 
the energy density and pressure of the partonic plasma, and c s is the sound velocity. The 
vacuum energy density e vac (bag constant) is introduced to account for the fact that the 
chromofield and partonic plasma can exist only in the perturbative vacuum. In numerical 
calculations e vac is fixed to the value 0.4 GeV/fm 3 . It is assumed that the chromofield energy 
density e/(r) and the plasma energy density e p (r) are functions of the proper time r only, 
defined in the interval T] t {r) < rj < T] p (t). The plasma energy density, e p (r) is found from 
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the hydrodynamical equation with a source term due to the field decay: 

^+(1 + c 2)£p = _^/ . (7) 
dr v ' r dr 

We integrate the Eqs. (4)-(6) until the time when the total pressure in the region between 
the slabs vanishes, i.e. B(t) = 0. After this time the slabs move with constant velocity. 
Obviously, the solution of equation (j4|) is 

T 

P a = Pa(r )^T- [ B(r)rdT, (8) 
r t J 

TO 

and the slab rapidity Y a (r) can be easily found from the definition P a = M a sinh (Y a — r] a ). 
Then, the slab trajectory, z a (r), is obtained from the relation z a = rsinh^ a . 

Below we present results for Au+Au and d+Au collisions at maximum RHIC energy 
y / SA/ r /v=200 GeV per NN pair (yn = 5.4). The nucleon-nucleon inelastic cross section was 



taken to be a nn =42 mb 16J. The baryon mean transverse momentum (p±), which 



determines the internal energy of baryonic slabs, was fixed at the value 1.0 GeV/c, in 
accordance with BRAHMS data [(3] for midrapidity region. 

We have considered several functional forms for the time dependence of the chromofield, 
resulting in different plasma production rates and baryon deceleration patterns. Here we 

r i ~ 4 

present results for the power law, e/(r) = e/(r ) 1 + 2 -^ 3i characterized by the decay 
time Td = 0.6 fm/c, as expected for the Schwinger-like decay mechanism [3(). This function 
is shown in Fig. 1 together with the time dependence of the plasma energy density as 
predicted by the hydrodynamical equation ([7]). According to our picture, the Quark-Gluon 
Plasma(QGP) is produced by the continuous transformation of the field energy into the 
quark- ant iquark and gluon pairs. Because of the delayed production, the plasma energy 
density is always smaller than the initial energy density of the chromofield. In accordance 
with previous calculations [3], we find that it reaches only about 22% of the latter for the 
Schwinger-like decay law. 

Figs. 2 shows how the time evolution of the slab rapidities Y Ptt (r) depends on the initial 
energy density of the chromofield, as well as on the baryon numbers of colliding slabs. The 
ideal equation of state (c^=l/3) was used for the quark-gluon plasma. We present results 
for several values of the parameter e (see eq. ((3])) between 0.5 and 2.0 GeV/fm 3 . Upper 
panels show results for the collision of two equal slabs with baryon numbers N p = N t =5.8, 
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FIG. 1: Evolution of the chromofield energy density (thick solid line) and QGP energy density 
(thin solid line) in units of the initial chromofield energy density ej(ro). Results are shown for the 
Schwinger-like (power-law) decay mechanism with = 0.6 fm/c. 

representing an average pair of slabs in a central Au+Au collision. The initial energy 
densities of the chromofield, e/(ro), range in this case from 16.8 (lower curves) to 67.3 (upper 
curves GeV/fm 3 . According to the BRAHMS data the mean baryon rapidity loss for the 
most central Au-Au collisions is (Sy) ~ 2.0. From the figures one can see that slab rapidities 
Yp it (r) drop rapidly at very early stages of the deceleration process, when the field is still 
strong. At later times, not only the decreasing chromofield (see Fig. 1) but also the plasma 
counter-pressure (left panels) cause an early termination of the deceleration process. As the 
result, the asymptotic slab rapidities have a low sensitivity to the initial chromofield energy 
density. Nevertheless, the observed rapidity loss can be achieved with the initial field energy 
density e/(T )=50-f- 70 GeV/fm 3 . On the other hand, it is likely that the plasma pressure 
in the vicinity of the slabs is not as strong as in the central region. Then the deceleration 
process will be less affected by the plasma, and resulting slabs' rapidities will be smaller. 
This is illustrated in the right panels showing calculations of the slab rapidities ignoring the 
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FIG. 2: Projectile (upper curves) and target (lower curves) slab rapidities as functions of proper 
time calculated for the power-law chromofield decay with r^=0.6 fm/c. Different pairs of curves 
correspond to different values of the parameter eo (indicated in the figure). Results are shown for 
two cases: a,b) equal slabs with N p = Nt = 5.8 representing a central Au+Au collision, and c,d) 
two different slabs with N p =2, N t =8.8 representing a central d+Au collision. Left and right panels 
show the calculations with and without the back reaction of the produced plasma, respectively. 

plasma pressure. In this case the required rapidity loss can be achieved with a much weaker 
chromofield, e/(r ) 30 4- 40 GeV/fm 3 . It is clear that the realistic situation will be in 
between of these two extremes. 

Low panels show results for an asymmetric collision of N p —2 and N t =8.8 slabs, which 
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can be interpreted as a deutron colliding with a center of a gold nucleus. The parameter eo 
is the same as before but now it corresponds to c/(to) between 8.9 and 35.5 GeV/fm 3 . One 
can see that the rapidity shifts are now different for the light and heavy slabs. Moreover, the 
rapidity lost by the smaller slab is significantly larger as compared with the bigger one. This 
is of course a direct consequence of the fact that equal forces generate a larger deceleration 
for a smaller body. It is interesting to note that the mean rapidity loss of a deutron-like 
nucleus in a d+Au collision is predicted to be about two times larger than the mean rapidity 
loss of the gold nuclei in a central Au+Au collision. 

In Fig. 3 we present the net-baryon rapidity distributions calculated for the central 
Au+Au collisions at maximum RHIC energy y / J/v r /v=200 GeV per iViV-pair. The calcu- 
lations were done by applying Gaussian weights for different color charges generated on the 
slabs [2]. Since the energy density of the chromofield is proportional to the square of the 
areal charge density, like in a capacitor, these weights are exponential in terms of the field 
energy density, i.e. P{ef) oc exp — , where (e/) = ej(r ) is the mean energy density of 
the field as parameterized in eq. ([3]). The introduction of the fluctuating field leads to a very 
interesting effect: although the zero field has a maximum probability, quite large values of 
the field, comparable with (e^), are possible with a certain probability too. The larger fields 
lead to larger decelerations of the slabs. As a result, the dNs/dy distribution has a long tail 
towards central rapidities. This is exactly what is observed by the BRAHNS collaboration 

The results shown in Fig. 3 were obtain without back reaction of the partonic plasma on 
the baryon deceleration. As one can see, in this case we can get a perfect fit of the data 
by choosing e ~ 1 GeV/fm 3 , which, according to eq. d3J), corresponds to the actual field 
energy density of about 35 GeV/fm 3 . Inclusion of the back reaction of plasma would require 
increasing of the initial chromofield energy density or /and the field decay time. 

On the basis of our model we can make predictions for the net-baryon dynamics for the 
future LHC experiments. We use eq. ([3]) to extrapolate the initial field energy density from 
the RHIC (^=200 GeV) to the LHC (^=5500 GeV) energy domain. Then we get the 
mean net-baryon rapidity loss between 2.7 and 5.5 units for the calculation with and without 
the back reaction of partonic plasma, respectively. 

In conclusion, the collective deceleration of valence (net) baryons was studied assuming 
that strong longitudinal chromofields are formed at early stages of an ultrarelativistic heavy- 
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FIG. 3: Net-baryon rapidity distribution in central Au+Au collisions (-/V part =380) at ^/sNN = 200 
GeV per AA-pair calculated for 3 different values of the parameter eo (indicated in the figure). 
Back reaction of the partonic plasma is not included in the calculation of the slab rapidities. Dots 
are experimental data of BRAHMS collaboration [6(. 

ion collision. We have solved classical equations of motion for baryonic slabs under the action 
of a time- dependent chromoelectric field. A simpl;e hydrodynamical consideration shows that 
the energy density of produced plasma, due to the delayed production, is always significantly 
lower than the initial energy density of the chromofield. It has been demonstrated that 
the net-baryon rapidity loss (8y) ~ 2 can be achieved with the initial energy density of 
the chromofield in the range of 40 to 70 GeV/fm 3 . The calculated net baryon rapidity 
distributions for central Au+Au collisions at RHIC energies are in good agreement with 
BRAHMS data. The mean net-baryon rapidity loss for central Pb+Pb collisions at LHC 
energies is predicted between 2.7 and 5.5 units. 

More detailed results, including the rapidity densities of net baryons and partonic plasma, 
calculated for different field decay patterns and centrality classes, will be presented in the 
forthcoming publication [ig| ] . 
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